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We present the results of a combined experimental and numerical study of the generation 
of internal waves using the novel internal wave generator design of Gostiaux et al. (2007). 
This mechanism, which involves a tunable source comprised of oscillating plates, has so 
far been used for a few fundamental studies of internal waves, but its full potential 
has yet to be realized. Our studies reveal that this approach is capable of producing 
a wide variety of two-dimensional wave fields, including plane waves, wave beams and 
discrete vertical modes in finite-depth stratifications. The effects of discretization by a 
finite number of plates, forcing amplitude and angle of propagation are investigated, and 
it is found that the method is remarkably efficient at generating a complete wave field 
despite forcing only one velocity component in a controllable manner. We furthermore 
find that the nature of the radiated wave field is well predicted using Fourier transforms 
of the spatial structure of the wave generator. 


1. Introduction 

The study of internal waves continues to generate great interest due to the evolving ap¬ 
preciation of their role in many geophysical systems. In t he ocean, internal wave s play an 
important role in dissipating barotropic tidal energy (see lGarrett fc Kunzell2007l for a re¬ 
view) , whereas atmospheric internal waves ar e an important means of momentum trans¬ 
port ( Alexander. Richter h Sutherland 20061) . In both the o cean and the atmosphere, 
intern al wave activity also impacts modern day technology (lOsborne. Burch fc Scarlet 
I978I) . Many unanswered questions remain, however, particularly regarding the fate of 
internal waves. For example: How much mixing do internal waves generate in the ocean, 
and via what processes? And at what altitudes do atmospheric internal waves break and 
deposit their momentum? The ability to reliably model internal wave dynamics is key to 
tackling important questions such as these. 

Internal waves come in a wide variety of forms. The simplest, a plan e wave, is the 
basis of many theoretical stu dies that provide fundamental insight (see iThornd Il987 . 
I998I: Dauxois fc Younel 19991 for instance), especially as any linear wave structure can 
be decomposed into independent plane waves via Fourier transforms. While plane wave 
solutions are the focus of many theoretical studies, laboratory experiments and field ob¬ 
servations reveal that in ternal waves generated by a localized sou rce, such as the tidal 
flow past an ocean ridge (Bell||l975: Martin. Rudnic k fc Pinkell 200(il) o r deep tropical con¬ 
vection in the atmosphere ( Walterscheid. Schubert fc BrinkmanI 2001 ). produce coherent 
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wave beams that radiate away from the generation site. In vertically-finite domains, 
such as the ocean, the internal wave field can be conveniently described using vertical 
modes, i.e. horizontally propagating and vertically sta nding waves, whose spa tial form 
is dictated by the vertical stratification, as discussed in Echeverri_e£a^ J^OO^ . Finally, 
localized solitary wave structures are also ubiquitous ( New fc Da Silvaf 2002h . 

Investigation of these different internal wave forms in laboratory exper iments has 
played a key role in internal wave research, starting with the pioneering work of iMowbrav fc Rarity 
I I967ll on the wave beams generated by an osc illating cylinder. S i nce th en, internal waves 
have been generated using a variety of means. [Pelisi fc Orlanski ( 1975ll performed an ex¬ 
perimental study of the reflection of nominally plane waves, produc ed by a paddle mecha¬ 
nism, from a density jump. A similar paddle mechanism was used bv llvev. Winters fc De Silva 
ll2000ll to study th e dissipation caused by internal wave breaking at a sloping boundary. 
iMaas et a,l\ ( 1997t l used vertical oscillations of a tank filled with salt-stratified water to 
parametrically excite internal waves, wh ich eventually focused onto inter nal wave attrac¬ 
tors in the tank. With the ocean in mind, Gostiaux fc DauxoisI ( 2007ll and lEcheverri et~^. 

I 2009ll produced internal waves by side-to-side oscillation of topography. 


The aforementioned experimental methods of internal wave generation have three in¬ 
herent shortcomings. Firstly, they produce wave fields that are invariant in one horizontal 
direction, and thus nominally two-dimensional. In this paper, we too restrict ourselves 
to the study of such situations, using z to refer to the vertical direction, antiparal¬ 
lel to the gravity field g = —ge^, and x to the horizontal direction; the possibility of 
generating three-dimensional wave fields using the novel generator is raised at the end 
of the paper. The second shortcoming is that, with the exception of towed topography 
( Baines fc Hoinkalll985 : Aguilar fc Sutherland!2006l l. ]3re-existing methods radiate waves 
in multiple directions rather than in a single direction. This is due to the dispersion re¬ 
lation for internal gravity waves, 

Lu^ = siix^ 9 , ( 1 . 1 ) 

which relates the forcing frequency, w, to the local angle of energy propagation with 
respect to the horizontal, 9, via the Brunt-Vaisala frequency, N = y^—gdzp/p, where 
p is the background fluid density. Since waves propagating at angles ±9 and n ± 9 
all satisfy dm for a given frequency ratio w/iV, a two-dimensional localized source, 
such as a vertically oscillating cylinder, generates internal waves propagating in four 
different directions. Propagation in two of the four directions can be suppressed by either 
providing oscillations along only one of the directions of p ropagation (i.e. 9 and n + 9) 

I Gavrilov fc Ermanvuk 1996t Ermanvuk fc Gavrilov 2008 1 or by using a paddle system 
at a boundary ( Pelisi fc Orlanskil 1975 1. These arrangements nevertheless still produce 
an undesirable second set of waves that must somehow be dealt with in an experiment. 
The third, and perhaps the most significant, shortcoming is that all pre-existing methods 
provide very limited, if any, control of the spatial structure of an internal wave field. 

A major advance in internal wave generation recently o ccurred with the design of a 
novel type of internal wave generator (IGostiaux et al.l 1200711 . This design uses a series of 
stacked, offset plates on a camshaft to simultaneously shape the spatial structure of an 
experimental internal wave field and enforce wave propagation in a single direction, as 
illustrated in figure [TJ The maximum horizontal displacement of each plate is set by the 
eccentricity of the corresponding cam, and the spatio-temporal evolution is defined by 
the phase progression from one cam to another and the rotation speed of the camshaft. 
So far, this novel configurati on has been used to study pl ane wave reflection fr om sloping 
boundaries ( Gostiaip] 120061 1. diffraction through a slit (Mercier et al. 2008h a nd wave 
beam propagation through nonuniform stratifications ( Mathur fc Peacockll2009ll . Despite 
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Figure 1. A schematic showing the basic configuration of a novel internal wave generator. 
Plates are vertically stacked on an eccentric camshaft. See text in section 2 for the definitions of 
the different lengths W, £ and A. The circular arrow at the top of the generator illustrates the 
direction of rotation of the camshaft, the thick vertical arrows show the corresponding motion of 
the wave form of the plates, and the dashed oblique arrows indicate the resulting local velocity 
field, and Vg indicate the direction of phase and group velocity, respectively. 


these early successes, however, there has been no dedicated study of the ability of this 
arrangement to generate qualitatively different forms of internal wave fields, and several 
important questions remain. For example, how does a stratified fluid that supports two- 
dimensional waves respond to controlled forcing in only one direction (i.e. parallel to the 
motion of the plates)? 

In this paper, we present the results of a comprehensive study of two-dimensional 
wave fields produced by different configurations of novel internal wave generators, and 
reveal that this approach can accurately produce plane waves, wave beams and discrete 
vertical modes. The results of experiments are compared with predictions based on the 
Fourier transforms of the spatial structure of the wave generator, which proves to be a 
very useful and simple tool for predicting wave fields, and numerical simulations, which 
allow investigation of the boundary conditions imposed by the generator. The material 
is organized as follows. Section [2] presents the experimental and numerical methods used 
throughout the study. The generation of plane waves is addressed in § [3l followed by the 
generation of self-similar wave beams and vertical modes in § |4] and § El respectively. 
Our conclusions, and suggestions for future applications of the generator, are presented 
in §[6l 


2. Methods 

2.1. Experiments 

Throughout this paper, we consider the case of wave fields excited by a vertically stand¬ 
ing generator with horizontally moving plates of thickness i, as depicted in figure [T] 
This scenario, which is possible because the direction of wave propagation is set by the 
dispersion relation CID, has two major advant ages over the other possibility of a genera¬ 
tor ti lted in the direction of wave propagation (jGostiaux et ad 120071 : iMathur fc Peacock 
2Q0^. First, it is far more convenient because it requires no mechanical components to 
orient the camshaft axis and no change of orientation for different propagation angles. 
Second, unwanted wave beams that are inevitably produced by free corners within the 
body of a stratified fluid are eliminated because the generator extends over the entire 
working height of the fluid. 
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Generator Height Width Plate thick. Plate gap Max. eccentricity Num. of plates 

Lyon 390.0 140.0 6.0 0.7 10.0 60 

MIT 534.0 300.0 6.3 0.225 35.0 82 

Table 1. Details of the wave generators at ENS de Lyon and MIT. Dimensions are in mm. 


Two differ ent laboratory experimental facilities, both using the double-bucket method 
I 0steiill965ll to create salt/density stratifications, were used. The first, at ENS de Lyon, 
utilized a 0.8 m long, 0.170 m wide and 0.425 m deep wave tank. The wave generator, 
whose characteristics are listed in table [TJ was positioned at one end of the tank. On 
each side of the wave generator there was a 0.015 m gap between the moving plates 
and the side wall of the wave tank. Visualizations and quantitative measurements of 
the densit y-gradient perturba tion field were performe d using the Synthet i c Sch lieren 
technique ( Dalziel et al. 2000f) . The CIV algorithm of Fincham fc Delei^ ( 2000f) was 
used to compute the cross-correlation between the real-time and the t = 0 background 
images. Blocksom filter matting was used to effectively damp end wall reflections of 
internal waves. The ENS Lyon setup w as used to run experimental s tudies of the classical 
Thomas-Stevenson wave beam profile ( Thomas fc Stevenson 1972ll . detailed in §|31 

The second system, at MIT, utilized a 5.5 m long, 0.5 m wide and 0.6 m deep wave 
tank. A partition divided almost the entire length of the tank into 0.35 m and 0.15 m 
wide sections, the experiments being performed in the wider section. The wave generator, 
whose characteristics are given in table [1] was mounted in the 0.35 m wide section of the 
tank with a gap of 0.025 m between the moving plates and either side wall. Parabolic 
end walls at the ends of the wave tank reflected the wave field produced by the generator 
into the 0.15 m wide section of the tank, where it was dissipated by Blocksom filter 
matting. Visualizations and quantitative measurements of the velocity field in the vertical 
midplane of the generator were obtained using a LaVision Particle Image Velocimetry 
(PIV) system. This facility was used for studies of plane waves and vertical modes, 
detailed in § |3] and § [5l respectively. 

Examples of the amplitude and phase arrangements of the plates for the experiments 
discussed in this paper are presented in figure [21 We use the following terminology: M 
is the number of plates per period used to represent a periodic wave form of vertical 
wavelength A, VE is the total height of the active region of the generator with nonzero 
forcing amplitude, A(z) is the eccentricity of a cam located at height z, and (j){z) is the 
phase of a cam set by the initial rotational orientation relative to the mid-depth cam 
{(j) = 0). The actual profile of the generator is given by Re , where Re stands 

for the real part. For a plane wave, A(z) is constant and (j){z) varies linearly over the 
active region of the generator. For the Thomas & Stevenson profile there is a nontrivial 
spatial variation in both A{z) and (f>iz) over the active region, and elsewhere A{z) is zero. 
Finally, for a mode-1 wave field, A{z) varies as the magnitude of a cosine over the entire 
fluid depth, while 4){z) jumps by tt at mid-depth. 


2.2. Numerics 

Complementary two-dimensional numerical simulations, in which excitation by the gen¬ 
erator was modeled by imposing spatio-temporal variations of the velocity and buoyancy 
fields along one boundary of the numerical domain, were performed. The simulations, 
which assumed a Newtonian fluid in the Boussinesq approximation, solved the incom- 
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(a) 


(b) 


(c) 



position of the cams for different profiles used throughout the paper. These include: plane waves 
for M = 12 and W = 2\ (o), Thomas & Stevenson beam (*) and a mode-1 internal tide (<). Thin 


lines drawn through the discrete points are the corresponding analytical forms being modeled. 


pressible continuity, Navier-Stokes and energy equations 

V • = 0, 


dtv + CV X v) X V = — Vq -b bez 




(2.1a) 

(2.16) 


and 


dfb -b (i? • V) 6 = —N^v • Bz -b (2-lc) 

where v = (u, w) is the velocity field, v the corresponding velocity magnitude, p = q — 
v^/2 the pressure, b the buoyancy field, related to the density by p = po (l — g~^N^z — g~^ 

where po is the density at z = 0, v the kinematic viscosity and k the diff usivity. _ 

The code used was an extension of that developed for channel flow (jGilbertlIl988 b 
to which the integration of the energy equation and the possibility of spatially varying 
an d time evolving bounda ry conditions on the plates were added, as already presented 
in Martinand et al\ ( 2006ll for thermal convection. The method proceeds as follows. A 
numerical solution in a rectangular domain [0, lx\ x [0, lz\ is obtained by a tau-collocation 
pseudo-spectral method in space, using Fourier modes in the z-direction and Chebyshev 
polynomials in the confined cc-direction. This method very precisely accounts for the 
dissipative terms. The nonlinear and diffusion terms are discretized in time by Adams- 
Bashforth and Crank-Nicolson schemes, respectively, resulting in second order accuracy 
in time. As the simulation focuses on the linear and weakly nonlinear dynamics, the 
de-aliasing of the nonlinear term in the spatial expansion of the solution is not of crucial 
importance and is omitted to decrease the computational cost. Owi ng to the assumption 
of a divergen ce-free flow, the influence matrix method, introduced in Kleiser fc SchumannI 
( 198(1 Il984ll . is used to evaluate the pressure and the velocity field; the pressure gradient 
in the Navier-Stokes equation then being discretized by an implicit Euler scheme. Finally, 
the buoyancy term is also discretized in time by an implicit Euler scheme, since the energy 
equation is solved before the Navier-Stokes equations. 

Simulations were run by imposing forced boundary conditions on components of the ve- 
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locity and buoyancy fields at x = 0; no forcing was applied to the pressure, since its value 
on the boundaries is an outcome of the numerical method. The governing equations m 
were thus integrated together with Dirichlet boundary conditions 

V (0, 0, t) = Vf {z, t), b (0, z, t) = bf {z, t) ; (2.2) 

while 

V {G,z,t) = 0, b{l^,z,t) = 0, (2.3) 

were applied dX x = We note that this numerical forcing is Eulerian in nature, whereas 
the corresponding experimental forcing is Lagrangian in spirit. The spectral method 
introduces periodic conditions in the ^-direction, which have to be accounted for to 
avoid Gibbs oscillations. Therefore, the boundary conditions (12.21) were multiplied by a 

polynomial “hat” function H (z) = ^1 — (2z/lz — 1)^°^ , vanishing at z = 0 and z = Iz- 

The choice of the exponents in H (z) is qualitative, the aim being that the variation of 
the profile envelope be smooth compared to the spatial resolution, yet sharp enough to 
keep a well-defined width of forcing. 

The boundary conditions (12.31) imply wave reflection, with the reflected waves eventu¬ 
ally interfering with the forced waves. Thus, the numerical domain was made sufficiently 
large to establish the time-periodic forced wave field near the generation location long 
before reflections became an issue. For a typical simulation, the domain was G = 3.01 m 
long and G = 1.505 m high, and the number of grid points used was Nx = 1024 and 
Nz = 512, giving a spatial vertical resolution of 2.9 mm that ensured at least two grid 
points per plate. Satisfactory spectral convergence was confirmed for this spatial resolu¬ 
tion and the time step was set to ensure stability of the numerical scheme. 


2.3. Analysis 

A detailed study of the i mpact of sidewall boundar y conditions on the generation of shear 
waves was performed bv lMcEwan fc BainesI ( 1974ll . Here, we take a simpler approach and 
show that a useful tool for investigating both theoretical and experimental internal wave 
fields produced by the novel generator is Fourier analysis. This allows one to decompose 
internal wave fields into constituent plane waves, and readily make predictions about the 
radiated wave field. 

For an unconfined, inviscid, two-dimensional system, any physical field variable asso¬ 
ciated wit h a periodic internal wave field of frequency ui can be described by its Fourier 
spectrum ( Tabaei fc Akvlas 2003 : Tabaei. Akvlas fc LambllioO^ . i.e. 


ip(x,z,t) = 




2n 


+ 00 


Q^{kx,kz)P^'^^^+'^^^U {{kl + kl)u;^ - N^kl) dkzdkx, (2.4) 


/ —<X)J —oo 


where 0 represents a field variable (e.g. 6, u) and the Dirac ^-function ensures the 

dispersion relation (EB is satisfied by the plane waves components. Propagating waves 
in a single direction, say towards positive x-component and negative z-component for 
the energy propagation, require 

Q^{kx.kz)=d Vfc, <0 & Vfc, <0, (2.5) 

as noted bv iMercier et al\ ()2008ll . 

At a fixed horizontal location xo, the values of 'ip(xo, z, t) for all the field variables can 
be considered as boundary conditions that force the propagating wave field ^/>(x, z,t). 
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Knowing the Fourier transform of the boundary forcing, 

/• + 00 


Qip ( 3^0 ; ) — 


\/^ J- 


I 

J —( 


il){xo,z,t)e 


( 2 . 6 ) 


leads to complete description of the radiated wave field for x > xq'- 

p—iojt ^+ 00^+00 

Ip{x,z,t) = - 


271 


QA^o,h) {{kl + kl)oo^ - N^kl) dfc, dfc,, (2.7) 


assuming that only right-propagating waves (i.e. kx > 0 ) are possible. 

In practice, the novel wave generator we consider forces only the horizontal velocity 
field in a controlled manner, i.e. 


'tp{0,z,t) = u{0,z,t) = Re{U{z)e . 


( 2 . 8 ) 


As such, we expect the Fourier transform of this boundary condition to act only as a guide 
for the nature of the radiated wave field, since it is not clear how the fluid will respond to 
forcing of a single field variable. Throughout the paper, we perform the Fourier transform 
along a specific direction using the Fast Fourier Transform algorithm. To compare spectra 
from theoretical, numerical and experimental profiles with the same resolution, a cubic 
interpolation (in space) of the experimental wave field is used if needed. 

Unless otherwise stated, the experimental and numerical results presented are filtered 
in time at the forcing frequency uj. The aim is to consider harmonic (in time) internal 
waves for which we can define the Fourier decomposition in (12.41) . and to improve the 
signal-to-noise ratio, which lies in the range 10 ^ — 10 ^, with the best results for shallow 
beam angles and small amplitude forcing. The time window At used for the filtering is 
such that wAt/27r > 9 for Aq = 0.005 m and a;At/27r > 4 for Aq = 0.035 m (where 
Aq is the amplitude of motion of the plates defined in l3.1|) . ensuring sufficient resolution 
in Fourier space for selective filtering. The recording was initiated at time tg after the 
start-up of th e generator such that Ato/27r ~ 30 ^ 1, ensuring no transients remained. 
(IVoi.sinll2003h . 


3. Plane waves 

Since many theoretical results for inter nal waves are obtained for plane waves (e.g. 
ThorDel[i987 . 1998HDauxois fc YounsJll9^ . the ability to generate a good approximation 
of a plane wave in a laboratory setting is important to enable corresponding experimental 
investigations. In order to generate a nominally plane wave, however, one must consider 
the impact of the different physical constraints of the wave generator, which include: 
the controlled forcing of only one velocity component by the moving plates, the finite 
spatial extent of forcing, the discretization of forcing by a finite number of plates, the 
amplitude of forcing, and the direction of wave propagation with respect to the camshaft 
axis of the generator. In this section, we present the results of a systematic study of the 
consequences of these constraints; a summary of the experiments is presented in table 


3.1. Analysis 

Two-dimensional, planar internal waves take the form ■0 {x, z, t) = Re |• 0 pe(l'=xa:-l-l/c 2 z-la;t) | 
where k = (kxykz) is the wave vector and ip{x,zA) represents a field variable. Being of 
infinite extent is an idealization that is never realizable in an experiment. To investigate 
the consequences of an internal wave generator being of finite extent, the horizontal ve¬ 
locity boundary conditions used to produce a downward, right-propagating, nominally 
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plane wave can be written as: 

m(0, z, t) = Re {[0(z + iy/2) - 0(z - Vb/2)] AqU} } , (3.1) 

where fee > 0 is the desired vertical wave number, Ae = 2njke the corresponding vertical 
wavelength, —W/2 < z < W/2 the vertical domain over which forcing is applied, Aq 
the amplitude of motion of the plates, and 0 the Heaviside function. The spatial Fourier 
transform of m is: 

Qu{0, kz) = How-^sinc , (3-2) 

with sinc(a:) = sinx/x the sine cardinal function. In the limit VF —>■ oo, Eq. (lO) ap¬ 
proaches a delta function, which is the Fourier transform of a plane wave. Owing to the 
finite value of W, however, <5^(0, k^) does not vanish for negative values of suggest¬ 
ing that will also excite upward propagating plane waves. Following the convention 
usual in optics that <3^(0, fcz) is negligible for — ke\W/2 > 7t, if kg » 2n/W then (12. Sp 
is reasonably satisfied. 

Another consideration is that the forcing provided by the wave generator is not spa¬ 
tially continuous, but discretized by Np oscillating plates of width i. Accounting for this, 
the boundary forcing can be written as: 

r jvp-i I 

u(0, z, <) = Re - fo) - - fo + ^)] \ , (3.3) 


where Zj = j£ — W/2 and Npi = W. The Fourier transform of (13.3p for the specific case 
of constant amplitudes A{zj) = Agw, Vj G {0, • • • , Ap — 1}, reduces to 


t I k £ 

Qu{0, kz) = Aow^=sinc ( 

V27t \ 2 


sm 


(fcz - ke)W 


sm 


(fcz - ke)£ 


(3.4) 


a classical result often encountered for diffraction gratings. Consequently, the magnitudes 
of W and £ in comparison to the desired vertical wavelength Ae = 2n/ke characterize the 
spread of the Fourier spectrum and the potential for excitation of upward propagating 
waves. 

In the following sections, we quantify the downward emission of waves using the pa¬ 
rameter /3d, defined as: 

/ \Qu{xo,kz)\^ Akz 

Pd = ^ -, (3.5) 

/ \Qu{xo,kz)\'^ dkz 

J —oo 

which is essentially the ratio of the total kinetic energy of the downward-propagating 
waves to the total kinetic energy of the radiated wave field. 


3.2. Configuration 

The MIT facility was used for these plane wave experiments (see table [T]). A variety of 
different configurations were tested and these are summarized in table The fluid depth 
was H = 0.56 ± 0.015 m and the background stratification was N = 0.85 rad s“^ for all 
experiments. An example configuration of the cams (amplitude and phase evolution) is 
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Expt. 

Experiment/Simulation 

Forcing 

M 

IF/Ae 

Aq (mm) 

e (deg.) 

Common case 

1 

exp./sim. 

partial 

12 

3 

5.0 

15 

Forcing 

2 

sim. 

complete 

12 

3 

5.0 

45 


3 

exp./sim. 

partial 

12 

3 

5.0 

30 

Angle 

4 

exp./sim. 

partial 

12 

3 

5.0 

45 

5 

exp./sim. 

partial 

12 

3 

5.0 

60 


6 

sim. 

partial 

12 

3 

5.0 

75 


7 

exp./sim. 

partial 

12 

2 

5.0 

15 

Width 

8 

exp./sim. 

partial 

12 

1 

5.0 

15 

9 

exp./sim. 

partial 

12 

2 

5.0 

45 


10 

exp./sim. 

partial 

12 

1 

5.0 

45 

Discretization 

11 

12 

sim. 

exp./sim. 

partial 

partial 

OO 

4 

3 

3 

5.0 

5.0 

15 

15 


13 

exp./sim. 

partial 

12 

3 

35.0 

15 

Amplitude 

14 

exp. 

partial 

12 

3 

35.0 

30 


15 

exp. 

partial 

12 

3 

35.0 

45 


Table 2. Summary of experiments and numerical simulations. M is the number of plates used 
for one wavelength, W/Xe is the spatial extent of forcing expressed in terms of the dominant 
wavelength, Aq is the eccentricity of the cams, and 6 is the energy propagation angle. For 
complete forcing, u, w and b were forced at the boundary. 


presented in figure!^ Plane waves were produced by configuring the Np plates of the wave 
generator with an oscillation amplitude Aq = 0.005 m, with the exception of experiments 
13 to 15 for which Aq = 0.035 m. Results were obtained for different forcing frequencies 
corresponding to propagating angles of 15, 30, 45 and 60°. Visualization of the wave 
field was performed using PIV, for which it was possible to observe the wave field in a 
40 cm-wide horizontal domain over the entire depth of the tank, save for a 1 cm loss near 
the top and bottom boundaries due to unavoidable laser reflections. The corresponding 
numerical simulations were configured accordingly. 

3.3. Results 

3.3.1. Forcing 

The consequences of forcing only a single component of the velocity field, which we 
call partial forcing, in comparison to forcing both the velocity field components and the 
buoyancy field (assuming they are related by the inviscid linear wave equation), which 
we call complete forcing, were investigated first using the numerical simulations. Here, we 
present the results of simulations performed using a sinusoidal boundary wave form with 
W = 3Xe, M = 12 and Ae = 78.8 mm (experiments 2 and 4 of table [2|). The magnitude 
and frequency of the boundary condition for horizontal velocity were Aq = 5.0 mm and 
w = 0.601 rad s“^, the latter giving 9 = 45°. 

Figures [3Ka) and (b) present snapshots of the horizontal velocity fields u produced by 
partial and complete boundary forcing, respectively, and there is excellent qualitative 
agreement between the two. More quantitative comparisons are provided in figure [3](c), 
which presents data along the vertical cuts C indicated in figures |3Ka) and (b), located 
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(a) 


(b) 


u{xc,z) 

X 10'^ X 10"^ 


(c) 




XX z 

Figure 3. Comparison of the numerically obtained horizontal velocity field, u, for plane wave 
beams forced by (a) complete and (b) partial forcing (experiments 2 and 4 in table [2j . (c) 
Horizontal velocity along the cut C in (a) and (b), located at xc = 0.05 m for complete (—) and 
partial (-) forcing. All lengths are in m and all velocities in m s“^. 


at xc = 0.05 m; this location was chosen because it is close enough to the generator that 
viscous damping has only an order 1% effec t on the wave field (based on the linear viscous 
theory for plane waves by iLighthil] ( 1978ll ), yet is sufficiently far from the generator to 
allow the wave field to adapt to the boundary forcing. The cross beam profiles in figure 
EIc) show the amplitude of u is roughly 20% lower for the case of partial forcing compared 
to complete forcing, but otherwise their forms closely match. Spectral information reveals 
no other discernable difference between the two, and both cases give Pd > 0.99, compared 
to the theoretical prediction of Pd = 0.97, revealing that almost all the energy is being 
emitted downward. Qualitatively and quantitatively similar results to those presented in 
figures EKc) and (d) were obtained for different physical quantities (e.g. w and b) of the 
wave fields, and for the other configurations of the generator listed in table [2] 


3.3.2. Angle of emission 

As one might expect, the quality of the wave field was best for shallower propagation 
angles, for which the horizontal velocity is a more defining quantity, and degraded for 
steep propagation angles, where w becomes the dominant velocity component. This is 
demonstrated in a qualitative manner by figures |3Ka)-(d), which present snapshots of 
the experimental horizontal velocity fields for 0 = 15, 30, 45 and 60°, corresponding to 
experiments 1, 3, 4 and 5 in tabled respectively. Although the forcing velocity u = Aquj 
is strongest for 9 = 60°, the wave field resulting from this partial forcing is not as strong 
and coherent as those at lesser angles. 

Another demonstration of the consequences of partial forcing is given in figureEl which 
presents the efficiency of the wave generator as a function of the forcing frequency for 
experiments and numerics. The efficiency is defined as the magnitude of the horizontal 
velocity component u (figure [Sf a)), or the velocity in the direction of wave propagation u' 
(figure [5Kb)), averaged over the central 0.1 m of the cut C across the wave field, compared 
to the forced horizontal velocity Aqu. For complete forcing, one expects the ratio ujA qoj 
to be one, whereas a smaller value indicates a less efficient mechanism. The results show 
that both ujA quj and u'jA qoj decrease with increasing uj/N. A simple physical argument 
for the decay of the response with the propagation angle could be that the generator 
provides an initial amount of kinetic energy that is redistributed by the flow into both 
horizontal and vertical motions. If all the energy is appropriately distributed, we expect 
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X X X X 

Figure 4. Experimental horizontal velocity fields u for different forcing frequencies correspond¬ 
ing to angles of propagation (a) 9 = 15° , (b) 30° , (c) 45° and (d) 60° . These are experiments 
1, 3, 4 and 5 in tabled respectively. All lengths are in m and all velocities in m s“^. 


u{xc)/Aouj (a,) u'(xc)/A qw 



and tilted (<) wave generator for the mean amplitude and standard deviation of (a) u/Aquj and 
(b) u' jAoui. Results are plotted as a function of (ui/N)'^ for different frequencies corresponding 
to 6^ = 15, 30, 45, 60 and 75° (numerics only). The solid lines correspond to u'/{Aocu) = 1 and 
the dashed lines to u'/{A quj) = cos 9. 


u' = AqU), and thus u = AqujcosO = Aquj^I — represented by solid lines in 

both panels of figure [3 Below (w/iV)^ = 0.5 the experiments and numerics follow this 
trend quite closely, but then depart from it significantly for higher frequency ratios. Also 
plotted as dashed lines in figure [5] is the relation u' = Aqui cos 9, implying u'/A qoj = 
xA- — and ujA^iw = 1 — which does a reasonable job of capturing the 

trend of the results, especially at higher frequency ratios. This relation implies that the 
energy associated with the motion of the plates along the direction of propagation is 
primarily responsible for setting the strength of the wave field. 

Finally, we analyze the evolution the wave field with increasing angle by computing Pd- 
For these experiments, the Fourier analysis in section [3T] predicts Pd = 0.97, independent 
of the angle of emission. Although the numerical values concur with this prediction, with 
Pd > 0.99 for all angles, we obtain values of 0.99 for 0 = 15°, 0.98 for 9 = 30° and 
45°, and 0.93 for 9 = 60° for the experiments. One possible reason for this decrease 
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in efficiency is the finite-amplitude lateral displacement of the plates of the generator, 
which could partially block the propagation of steeper waves; this is not taken into 
account in the numerics. Another possible reason is that at higher frequency of forcing 
the Reynolds number for the oscillating plates is larger, increasing the likelihood of more 
complex dynamics near the oscillating plates and thereby weakening their coupling to 
wave generation. These issues are raised again in section 13.3.51 

Overall, the results in this section and the previous section reveal that partial (horizon¬ 
tal) forcing by a vertically standing generator works well for 0 < 45°. This gives the user 
the freedom to perform experiments over a range of angles without having to re-orient 
the generator. For larger angles, however, it would seem prudent to use a generator with 
its perpendicular axis tilted toward the direction of propagation. This was confirmed by 
a series of experiments similar to cases 1 and 3 to 5 in table [21 but with the generator 
tilted at 15° to the vertical; these results are also included in figure O For angles smaller 
than 45°, there is almost no difference in the efficiency of the generator, but we found 
that the 60° wave field produced by the tilted generator was noticeably stronger and 
more coherent than that produced by the vertically-standing generator. 

3.3.3. Finite extent 

The Fourier analysis in sectionjOlpredicts that a consequence of a generator inevitably 
being of finite vertical extent is the production of undesirable waves that propagate in 
the vertical direction opposite to the principal wave field. This can be seen in both 
the numerical and experimental wave fields in figures [3] and SI which contain a weak, 
upward-propagating wave field in addition to the principal downward-propagating wave 
field. Experiments were therefore performed for W/Xe = 1, 2 and 3 to investigate how 
the strength of the undesirable wave field was influenced by the vertical restriction of the 
forcing, and to determine how well this was predicted by simple Fourier analysis. Except 
for varying W/Xe and having 9 = 15° , the configuration was the same as in the previous 
subsection. The experiments are listed as experiments 1, 7 and 8 in table [21 

A direct comparison of experimental and numerical horizontal velocity fields for W/Xe = 
3 is presented in figure [HI demonstrating very good agreement between the two, provid¬ 
ing confirmation that our numerical approach can reliably model the horizontal forcing 
provided by the plates. The profiles presented in figure |6Kc), obtained at the vertical 
cuts indicated in figures |6Ka) andEKb), have only one small, but noticeable, difference: 
slightly higher peaks at either end of the experimental velocity profile. A similar level of 
agreement was obtained for the vertical velocity profile. 

The normalized experimental spatial Fourier spectra Qu{xc,kz) at xc = 0.05 m and 
the normalized theoretical spectra Qu (0, are presented in figure jT] A standout fea¬ 
ture of the results is that the theoretical Fourier transform does a remarkably good job 
of predicting the spectrum of the experimental wave field. For W/Xe = 1, the spectrum 
is broadly centered around the expected vertical wavenumber ke = 79.9 m“^, and this 
principal peak becomes increasingly sharp for W/Xe = 2 and 3. This evolution is quanti¬ 
fied by the half width 5k, defined as the width of the principal spectral peak at half peak 
amplitude, the values of which are listed in table |3| for the three different configurations. 

Another notable feature of the spectra is that strength of the upward propagating 
wavefield {kz < 0) significantly decreases with increasing W/Xe- Computing the param¬ 
eter fid for experimental and numerical cases 1, 7 and 8 of table |2] quantifies this trend. 
As seen in table |3l 98% of the energy propagates in the desired direction for W/Xe > 2. 
Furthermore, we computed 5k/ke and fid for cases 4, 9 and 10 of table |2] too, and the 
results presented in table |3| show that the influence of the width is more significant than 
the influence of the angle of emission. 
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u{xc,z) 



3 ^ X Z 

Figure 6. Horizontal velocity fields for (a) experiment and (b) simulation for experiment 1 in 

Tabled (c) Vertical profiles along the cut C located at xc = 0.05 m for the experiment (-) 

and the simulation (—). All lengths are in m and all velocities in m s“^. 


Qu {x, kz) / maxfc^ Qu (x, kz) 



Figure 7. Comparison between spatial spectra Qu { x , kz ) computed from the experiments (—) 
at station xc = 0.05 m and the theoretical expression ( —) computed from on the boundary 
X = 0, for IF = Ae (top), W = 2Ae (middle) and W = 3Ae (bottom). 
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Sk/k, (15°) 

Pd (15°) 

Sk/ke (45°) 

Pd (45° 

) 

W/Ae 

num. 

exp. 

theo. 

num. 

exp. 

theo. 

num. 

exp. 

theo. 

num. 

exp. 

theo. 

1 

1.67 

1.11 

1.21 

0.97 

0.94 

0.94 

1.31 

1.08 

1.21 

0.97 

0.92 

0.94 

2 

0.58 

0.57 

0.60 

0.99 

0.97 

0.96 

0.61 

0.58 

0.60 

0.99 

0.96 

0.96 

3 

0.37 

0.39 

0.40 

0.99 

0.99 

0.97 

0.40 

0.39 

0.40 

0.99 

0.98 

0.97 


Table 3. Relative half width, Sk/k^, and relative energy of the downward propagating wave, 
Pd, for cases 1, 7 and 8 of table [5] corresponding to a propagation angle of 15°, and experiments 
4, 9 and 10 for a propagation angle of 45°. 


3.3.4. Discretization 

To study the impact of spatially discretized, rather than continuous, forcing, experi¬ 
ments were performed for WjX^ = 3 with M = A and M = 12; these being experiments 
12 and 1 in table [2l For comparison, corresponding numerical simulations were also per¬ 
formed for these two configurations, as well as for the idealized case M —i oo, which is 
listed as experiment 11 in table [5] and corresponds to forcing discretized on the scale of 
the grid resolution in the numerical simulations. 

Snapshots of the experimental and numerical wave fields for M = A are presented in 
figuresIHl^a) and|8l[b), respectively, while figure|8l[c) presents vertical cuts of the horizontal 
velocity field at xq = 0.05 m for these two data sets. Even for this coarse discretization, 
there is still a remarkably smooth and periodic wave field that looks little different to that 
obtained using M=12 (see figure [i). And once again there is good agreement between 
experiment and numerics, with the slight exception of the outer edges of the profile where 
the numerical peaks are a little larger amplitude. 

Although the cross-section of the emitted downward-propagating, nominally plane 
wave looks reasonable, the discretization does induce more of an undesired upward- 
propagating wave, which can clearly be seen in figures ISJa) and jS^b) . Fourier spectra for 
experiments, numerics and theory corresponding to M=4, 12 and oo are presented in 
figure [9] and, as predicted by (13.4L the strength of the negative wave numbers noticeably 
increases with decreasing M. Most notably, for M = 4 (corresponding to £ = 19.6 mm) 
the discretization introduces a peak around kz = —235 (±4) m“^, which is strongest in the 
theoretical spectrum but nevertheless evident in the experimental and numerical spectra. 
By analogy with the theory of optical gratings, this value is in good agreement with the 
canonical formula 27r/Ae — 2tt/£ = —241 m“^, which can also be inferred from (13.41) when 
i <C Ae. For all three cases the principal peak remains sharp, with 5k/ke = 0.42. The 
value of Pd is 0.96 when M = 4, so a vast majority of the energy is still in the downward 
propagating wave field. 

3.3.5. Amplitude 

All the results presented thus far have been for Aq = 5.0 mm. To investigate the impact 
of a significantly larger amplitude of forcing on the quality of the radiated wave field, we 
performed a series of experiments with the same parameters as experiments 1, 3 and 4, 
with the exception of Aq = 35.0 mm; they are listed as experiments 13 to 15 in table [2j 

We found that the qualitative level of agreement between experiment and numerics 
for snapshots of the wave field was comparable to that presented in figure IHKa) and 
(b). When a more quantitative comparison is made, however, some consequences of 
the higher-amplitude forcing become apparent. For example, for vertical cuts located 
at xc = 0.075 m the amplitude of the horizontal velocity component in the experiments 
was 4.4 ± 0.9 mm s“^, compared to 5.50 ± 0.25 mm s“^ in the numerical simulations. 
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XX z __ 

Figure 8. Horizontal velocity fields for (a) experiment and (b) simulation for case 12 in Table[2] 

with M = 4. (c) Vertical profiles along the cut C located at xc = 0.05 m for the experiment 
(-) and the simulation (—). All lengths are in m and all velocities in m s“^. 


Qu{x, k^)/ maxfc^ Qu{x, k^) 



computed from (EU on the boundary a; = 0, for M = 4 (top), 12 (middle) and oo (bottom). 


We also note that although the forcing amplitude was increased by a factor of 7, in the 
experiments the wave amplitudes only increased by a factor of around 5. 0 

t We had to make these, and later, comparisons, for a vertical cross section further away from 
the generator than in our previous studies because the much larger amplitude motion of the 
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by their maximum amplitude. The data was obtained from cuts at xc = 0.075 m for (a) 
^0 = 0.005 m and (b) Aq = 0.035 m, for experiments (-) and numerics (—). 

Figure [TO] presents vertically-averaged temporal spectra of the horizontal velocity for 
the cuts at xq = 0.075 m, for both small and large amplitude forcing. The numerical 
data, being at somewhat higher temporal resolution, has a lower noise level than the 
experimental data. Both experimental and numerical spectra display the same qualitative 
change; the large amplitude forcing introduces more significant higher harmonic content 
into the wave held. 

The normalized spatial spectra of the wave held for frequencies corresponding to 9 = 
15, 30 and 45° are presented for both small and large amplitude forcing in hgure[TT] The 
angle of emission does not seem to signihcantly impact the quality of the wave held for 
the small amplitude forcing, but this is not so for the larger amplitude forcing, for which 
we hnd that jdd decreases from 0.97 to 0.94 for 9 = 15° and 30°, respectively (although 
5k/ke = 0.38 remains constant). For 9 = 45° and Aq = 35.0 mm, the generator no 
longer generates a clean plane wave, the main peak being centered around 41.5 m“^ 
and Pd = 0.80. The cause of this breakdown is not easy to discern. One hypothesis is 
that breakdown occurs at criticality, when the angle of wave propagation exceeds the 
maximum slope angle of the face of the generator, 0m = § — arctan ^ ^ . The reason 

is not as simple as this, however, since 0m = 20° for the experiments with Aq = 35.0 mm, 
and yet the generator is still an efficient source of plane waves for 0 = 30°. Other factors, 
such as the nonlinear coupling between the plates and the wave field, as characterized 
by the Reynolds number of the plate motion, would also seem to play a role. 

3.4. Summary 

Through a systematic series of experiments, listed in Tabled we can draw several conclu¬ 
sions about the ability of the novel wave generator to generate plane waves. We find that 
the spatial Fourier transform of the profile of the wave generator can reasonably predict a 

plates created a more intense wave field very close to the generator, where it was not possible 
to get reliable experimental data. 
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Qu{x, fcz)/maxfe^ Qu{x, kz) 



Figure 11. Comparison between spatial spectra Qu{x,kz) computed from the experiments 

with Aq = 0.035 m (x) and j4o = 0.005 m (-) at station xc = 0.075 m and the theoretical 

expression (—) computed from (13.41) on the boundary a; = 0, for 9 — 45° (top), 30° (middle) 
and 15° (bottom). 


priori the spectrum of the radiated wave field. If more comprehensive resources are avail¬ 
able, a numerical simulation with boundary forcing applied at a;=0 can reliably reproduce 
the emitted wave field for small-amplitude forcing. As one might expect, the spectrum 
of the wave field becomes increasingly sharp about the dominant wavelength, and thus 
more akin to a plane wave, as the number of wavelengths excited increases, and even a 
very crude spatial discretization of the desired wave form produces a remarkably smooth 
and coherent wave field. For large amplitude forcing, the main impact is an increase of 
the harmonic content of the wave field. Overall, we conclude that a vertically-standing 
wave generator produces a radiated wave field of high quality provided 0 < 45°. 


4. Wave Beams 

Wave be ams are a common feature of internal wave field s in both laboratory ex- 
periments (Peacock fc Tabei 2005 : Gostiaux fc Dauxois 20071) a nd geophysical settings 
( Lam. Maas fc Gerkema 2004 Martin. R.udnick fc Pinkell 2006[l . since they are read¬ 
ily generated by periodic flow relative to an obstacle, be it a cylinder or an ocean 
ridge for example. We cho ose to investigate the so-called Thomas-Stevenson profile 
( Thomas fc StevensonI 1972ll . a viscous self-similar solution of (EH) that can be con¬ 
sidered as the fa r-field limit of the vis cous solution of elliptic cylinder oscillating in 
a stratified fluid (|Hurlev fc Keadvl 119971) . It has been shown that this profile describes 
oceanograp hically relevant internal wav e beams far from their generation site at the conti¬ 
nental shelf (jGostiaux fc Dauxoisll2nn7n : after their reflection at the bottom of the ocean, 
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suc h wave beams are thought to be the c ause of solitons generated at the thermocline 
(see New fc Da Silva 2002 : Gerkema 2001 . for instance). 


4.1. Analysis 

Consider a downward-propagating, right-going beam at angle 9 G [0,7r/2] with respect 
to the horizontal, here 9 being defined to be positive in a clockwise sense. Let ^ = 
X cos 9 — zsin9 + I and rj = xsm9 + z cos 9 be the coordinates parallel and transverse to 
the wave beam, respectively, with I corresponding to the distance from the point source 
to the origin of the cartesian frame at the center of the active region of the generator. 
At leading order, the parallel and transverse velocity components and buoyancy fields of 
the Thomas-Stevenson profile are: 


u' = uo 


sin 9 


-2l3 

Re 


|y fcexp (—exp — iwt^ dfc| , (4.1a) 


v' (^,?7,t) = uo 


sin 9 
V .9 




- 2/3 

Re 


—i J exp {—k^) exp — iwtyfc| , 

(4.1&) 


and 


^ (Cj 9; 0 =-^^0 ^^ Re-|—i J fc exp (—exp — iwt^ dfel, 


where uq is the amplitude of the horizontal velocity and a = (2Ncos9/. 

In principle, to most accurately reproduce 63 ) one should tilt the generator and 
configure the profile of the forcing plates to match the transverse profile of the Thomas- 
Stevenson beam. As stated earlier, however, we consider a vertically-standing wave gen¬ 
erator, since one will typically want to investigate several wave beam angles in an ex¬ 
periment and reconfiguring the system for each angle is impractical. Thus we seek to 
reproduce the profile (j4.1l) by using only the longitudinal velocity profile (14.lap , instead 
of the true horizontal velocity profile, to prescribe the forcing at a;=0, i.e. 


u (0, 0 , t) = u {I, z, t). 


(4.2) 


This is an approximation of the exact solution, which will become increasingly valid with 
decreasing 9. For a given viscosity, stratification and forcing frequency that determine a, 
equation (14.21) sets effective values for the parameters ^ and mq in (14.11) . 


4.2. Configuration 

The experiments were performed using the ENS Lyon generator (see table [T]) with a 
background stratification N = 0.82 rad s“^ and forcing frequencies oj = 0.20 rad s“^, 
0.44 rad s“^ and 0.58 rad s“^, corresponding to propagation angles of 14, 32 and 44° 
respectively. The arrangement used twenty-one plates to discretize the profile at a: = 0, 
thirteen of which covered the 9 cm active region. Based on the results of section 13.3.41 
this level of discretization is expected to be sufficient to resolve the structure of the wave 
beam. The configuration of the cams (amplitude and evolution of the phase) is depicted 
in figured with the maximum amplitude of oscillation being 10 mm. 

Experimental visualizations were performed using the Synthetic Schlieren method, 
which gives direct measurements of the gradient of the buoyancy field. For the following 
study, we integrate this data and compare the measured buoyancy field with the analytical 
model (j4.1cl) and numerical results. For these simulations, the numerical domain was 
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0 0.05 0.1 0.15 0.2 0.25 0 0.05 0.1 0.15 0.2 0.25 


Figure 12. Buoyancy field b{x,z,t) for a 14° Thomas-Stevenson beam. Black lines are the 
transverse cuts used in figure[T3]to compare with the model, (a) Experimental and (b) numerical 
buoyancy fields are in rad.s”^. 


0.80 m long and 0.60 to 1.01 m high, with a vertical resolution Az = 0.79 mm. The 
forcing on the boundary was discretized on the scale of the numerical grid. 


4.3. Results 

A direct comparison between experimental and numerical buoyancy fields for 6 = 14° is 
presented in figure [TH There is close qualitative agreement between the two, and it is 
notable that there is no visible sign of any upward propagating beam coming from the 
generator, due to the highly resolved, smooth boundary conditions that were used. For 
a vertical cut at xc = 0.05 m of the vertical component of the density gradient, d^h, we 
find that Pd = 0.98 for the experiments and Pd = 0.99 for the numerics. We note that the 
maximum values of dzb are around 2% of the background stratification, and so the wave 
field can reasonably be considered linear. Analysis of temporal spectra confirmed that 
harmonics were at least one order of magnitude smaller than the fundemantal signal. 

For a more quantitative comparison with the self-similar solution of Thomas-Stevenson, 
we investigate the buoyancy field extracted along the cuts indicated in figures HHa) and 
(b). Specifically, we consider the normalized transverse profiles: 


bmiO 


poo 

/ k exp 1 

- 1 - ikar]/^^^^^ 

1 dk 

J 

poo 

0 

k exp {—k^) dfc 



(4.3) 


where 

^m(f) = niax 6 (^, ?7, 0) = A^uq ^ J kexp[—k^)dk (4.4) 

is the maximum amplitude of the buoyancy perturbation (I4.1c|) along a transverse cut, 
which lies at the center of the beam (77 = 0 ). The results are presented in figure [TSl' al 
and (b), where it can be seen that both the experimental and numerical results evolved 
spatially in a self-similar manner, with only small differences compared to the analytical 
model. 

We also confirmed the ability of the generator to produce waves beams for some steeper 
angles, by performing experiments for 6 = 32° and 6 = 44°. Comparisons with the 
analytical solution revealed the same level of agreement as for the 9 = 14° wave beam. 
The wave generator still emitted a beam in only one direction, and the experimental 















20 


Mercier, Martinand, Mathur, Gostiaux, Peacock and Dauxois 


(a) (b) 



experiments and (b) numerical simulations as a function of for the case 6 = 14°. Aster¬ 

isks correspond to the different cuts in figure m\ the dashed line corresponds to the analytical 
solution in 


values of (3d were 0.99(8) for 9 = 32° and 0.99(8) for 6 = 44°, with no discernable change 
in the value of 13d for the numerical simulations. 

4.4. Summary 

On the basis of this study, we conclude that the novel wave generator is capable of 
producing a wave beam structure of a desired form, which demands excitation of a 
prescribed Fourier spectrum. We have demonstrated this for the example of the Thomas- 
Stevenson profile, and speculate that although it is not necessary, perhaps even closer 
agreement with theory can be obtained using a tilted generator. 


5. Vertical modes 

Vertical internal wave modes play an important role in our current understanding of 
internal tides in the ocean (IGarrett fc Kunzell2007ll . To date, however, there has been lit¬ 
tle pro gress in producing high-quality vertical modes in laboratory experiments. iThoroel 
1 1968h generated a mode-1 disturbance b y oscillating a flap hinged about a horizontal axis 
at mid-depth; and iNicolaou et cH. ( 1993 1 produced low modes in an essentially two-layer 
system with a thermocline, but neither of these approaches can readily produce arbitrary 
modes in an arb itrary stratificat i on. G eneralized forcing of a spectrum of vertical modes 
was obtained bv lEcheverri et al ( 2009ll . who used an oscillating Gaussian topography to 
generate an internal wave field, and developed a robust algorithm for extracting modal 
amplitudes from experimental data. Here, we demonstrate the ability of the novel wave 
generator to reliably produce arbitrary internal wave modes. 

5.1. Analysis 

For a stratified fluid of constant TV, the horizontal velocity field associated with the n**' 
vertical mode of frequency oj propagating from left to right is: 


Un{x, z, t) = Re 


uo cos 


/riTTZN 


exp 1 1 


0 : 


nn 


• X — \(jjt 


(5.1) 


\ H J K'Hcote' 

where uq is a complex amplitude that sets both magnitude and phase, n is a positive 
integer, z = 0 and z = H are the bottom and top boundaries, respectively, and 9 is 
the first-quadrant angle that satisfies the dispersion relation (HH). Note that 9 does not 
specify the direction of energy prop agation fo r a vertical mode and only plays a part in 
setting the horizontal wavenumber ( Gilll 198^ . 
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The idealized boundary forcing of the horizontal velocity at a; = 0 to excite the n**' 
vertical mode is: 

(5.2) 


( TllXZ \ 

cos(wt), 


where a is an arbitrary amplitude. To analyze the quality of the wave field generated by 
this forcing, instead of Fourier transforms one uses modal analysis, which is equivalent 
to Fourier series in a constant stratification. The horizontal velocity component of the 
resulting wave field can be written as: 


u{x,z,t) = ^ 


n —1 


(nnz\ ( nnx 

an cos y —j cos 


— cut -|- ( 


.)■ 


(5.3) 


where a„ and (j)n are the strength and the phase of the mode respectively. Similar 
results exist for other physical variables, including vertical velocity and buoyancy fields. 

From a practical point of view, one decomposes the experimental generated wave field 
at a fixed x location into the vertical basis m odes using the numerical algorithm described 
and implemented in Echeverri et al. ( 20091) . The modal decomposition is then performed 
at several other x locations, and the variations in a„ and cpn across various x locations 
gives an estimate of the experimental errors in the results. One can reliably correct for 
viscous dissipation of the modes, if needs be, by introducing in (15.31) a multiplicative 
term of the type with 


(5.4) 


vn^ f 


2lo \h) 



being the spatial damping rate ( Thoroel 1968t Echeverri et al. 20091) . In our experiments, 
fn is very small (roughly 10“"'’ m“'^ for mode-1 and lO”"^ m“'^ for mode-2) and hence 
viscous dissipation can be neglected. 

The discretization results for plane waves in section 13.3.41 suggest that there could be 
limitations on the ability to resolve a vertical mode due to the discretization of the wave 
generator. By computing the modal decomposition of the discrete plate arrangement, we 
find that provided at least 3n equispaced plates are used to represent the mode, then 
more than 95% of the boundary forcing is contained in the mode. To ensure good 
quality of the generated wave field also requires one to account for the approximation 
that the boundary forcing of horizontal velocity occurs at a fixed x location even though 
the plates are actually moving. For this approximation to hold, the ratio of the maximum 
amplitude of oscillations of the plates to the horizontal wavelength corresponding to the 
n**' mode should be much smaller than unity, i.e. nA^/{2H coiO) <C 1. Provided these 
two conditions are satisfied, a high-quality wave field is to be expected. 


5.2. Configuration 

The MIT facility was used for the experiments with Np = 64 of the total 82 plates 
(see table [1]) . Individual vertical modes were produced by configuring the Np plates of 
the wave generator to reproduce (15.21) . The amplitude of oscillation of the plate 
centered at vertical position zj = jT/2 of mode n was A{zj) = Aq cos{nnzj/H), with £ 
being the plate thickness (see Fig. [2] for more details of the configuration). Experiments 
were performed for modes I and 2. The spatial resolution of the forcing, relative to the 
vertical wavelength of the mode being forced, was 1/64 in both experiments. The fluid 
depth was H = 0.416 m, the maximum amplitude of oscillation was Aq = 2.5 mm and 
the stratification was N ~ 0.85 rad s“'^. Visualization of the wave field was performed 
using PIV technique. Using this arrangement it was possible to observe the wave held in 
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Figure 14. Snapshots of experimental velocity fields for (a) mode-1 and (b) mode-2. The 
location of the wave generator is a: = 0 and the vertical dashed lines bound the domain over 
which modal decomposition is performed. The x and 2 coordinates are in meters and the grey 
scale is the velocity magnitude in m s“^. Arrows indicate local velocity direction. 


a 45 cm-wide horizontal domain and covering the full depth of the tank. No corresponding 
numerical simulations were performed in this case. 


5.3. Results 

We present detailed results for mode-1 and mode-2 wave fields with 6^=45°, for which 
the horizontal wavelengths excited were kx = 7.55 m“^ and kx = 15.1 m“^, respectively. 
Snapshots of the wave fields obtained in the vicinity of the generator are presented in 
figure [m One can clearly see the characteristic structure of a single vortex that covers 
the entire vertical domain for mode-1 ffigure rCT all. whereas for mode-2 the structure 
comprises stacked pairs of counter rotating vortices (figure [TU^b)). 

Modal decomposition of the experimental wave fields was performed at 81 ^-locations 
in the regions bounded by the vertical dashed lines in figures fTdl al and ITdlbl , and the 
results are presented in figure 1151 Since the maximum values of both u and w were the 
same for these 9 = 45° wave fields, reliable values of a„ and (j)n were obtained from 
both components of the velocity field. The two wave fields were clearly dominated by 
their mode number, with by far the largest detectable amplitude being for modes-1 and 
modes-2 in the two respective experiments, and with very little variability across the 
experimental domain, emphasizing the weak impact of viscosity in these experiments. 
We also observed very little variability in the phase of the dominant mode across the 
visualization window, impl ying that the wave field s were highly spatially coherent. Since 
energy flux scales as a^/n ( Echeverri et all ( 20091) 1 over 98.8 % of the energy was in the 
desired mode in each experiment. The efficiency of conversion, defined as anjAuj, was 0.89 
and 0.78 for the two experiments. The insets in figure [TS] present the vertically-averaged 
temporal Fourier spectrum of u. These are dominated by the fundamental frequency, 
demonstrating that very little higher-harmonic content was generated by nonlinearity. 


5.4. Summary 

Overall, these experiments reveal that the novel wave generator is capable of producing 
very high-quality radiating vertical modes in a linear stratification. For smaller angles 
(i.e. lower frequencies), 9 = 15° for example, we found that the wave field took longer 
horizontal distances to evolve into established modal solutions, presumably due to the 
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Figure 15. Modal decompositions of the wave fields presented in (a) figure fTlT a'). and (b) figure 
lldl bl. Error bars are the standard deviation of modal amplitude for the 81 vertical cross-sections 
studied between the dashed lines in figure [Td] The insets show vertically averaged absolute value 
of the temporal Fourier spectra of u at a; = 10 cm and x = 5 cm in (a) and (b), respectively. 

longer horizontal wave lengths of the modes. For higher frequencies, 0 = 60° for example, 
we observed high-quality modal solutions within our visualization window, but with a 
small variation in the dominant modal strength for mode-2, possibly due to nonlinear 
effects. For higher angles, therefore, it would seem prudent to decrease the maximum 
amplitude of oscillations and/or reduce the stratification N (allowing for lower forcing 
frequency) in an effort to produce as clean a wave field as possible. In principle, this 
approach can be extended to experiments with nonlinear density stratifications, pro¬ 
vided the stratification is known a priori so that the generator can be configured for the 
appropriate vertical structure of the horizontal velocity field. 


6. Conclusions 

Through combined experimental, numerical and theoretical studies we have demon¬ 
strated that the novel type of internal wave generator, comprising a series of stacked, 
offset plates, can reliably shape the spatial structure of an experimental internal wave 
field and enforce wave propagation. This approach is similar in spirit to multiple- paddle 
techniques that have been developed for generating surface waves (see Naitol 20061 for in¬ 
stance). We have demonstrated the ability of the generator to produce three qualitatively- 
different types of wave field: plane waves, wave beams and vertical modes. This new 
technology therefore provides a very useful tool to study all manner of internal wave sce¬ 
narios in the laboratory, in order to gain insight into geophysically important problems. 
Furthermore, our studies reveal that the Fourier transform of the spatial profile of the 
wave generator provides a reasonably accurate prediction of the form of emitted wave 
field, making it a useful tool when designing experiments. 

There are numerous examples of where this new found capability can now be uti- 
lized. For example, ob s ervations in several loc ations, in particular in the Bay of Biscay 
( New fc Pingreel 1991)1 : New fc Da Silval 2002l l. have determined that an internal tidal 
beam striking the thermocline is responsible for the generation of solitons. A full under¬ 
standing of the generation mechanism has yet to be achieved, however, and laboratory 
experiments using the Thomas-Stevenson beam profile impinging on a thermocline could 
provide significant insight. Indeed, the interaction of wa ve beams with nonlinear f eatures 
in the density stratification is of widespread interest ( Mathur fc Peacock 200^, since 
this is also relevant to how and where atmospheric internal waves break and deposit 
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their momentum ( Nault fc Sutherland! 2007 ). In regards to ocean mixing problems, an 
open issue is to determine th e fate of the internal tide which, among other thing s, can 
be scattered by topography ( Johnston &: Merrifieldl 2003t iRav fc MitchumI 1997 ). The 
ability to directly generate vertical modes provides a new capability to study these im¬ 
portant processes in controlled settings. Other interesting avenues for research are the 
generation of shear waves (McEwan & Baines 1974) and extensions to three-dimensional 
wave fields, which could be achieved by introducing some horizontal spatial structure (in 
the y-direction) to the leading edge of the moving plates. 
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